On the possibility of a long range proximity effect in a ferromagnetic nanoparticle. 
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We study the proximity effect in a ferromagnetic nanoparticle having a vortex magnetization 
pattern. We show that for axisymmetric system consisting of a circular particle and a magnetic 
vortex situated at the center of it no long range superconducting correlations are induced. It means 
that induced superconductivity is localized in the small area near the superconducting electrode. 
However, in the real systems axial symmetry can be broken by either a shift of the magnetic vortex 
from the origin or geometrical anisotropy of the ferromagnetic particle. In this case a long range 
proximity effect is possible. 

PACS numbers: 



I. INTRODUCTION 

Proximity effect in hybrid ferromag- 
netic/superconducting (FS) structures reveals a reach 
physics originating from the interplay between magnetic 
and superconducting types of ordering (see Reful for 
review) . There are two essential features of the proximity 
effect in FS structures which make it different from that 
in superconductor/ normal metal (SN) structures. In SN 
structures the penetration length of a condensate wave 
function into the normal metal is determined by the 
normal metal coherence length = y^D/2TTT, where D 
and T are the diffusion coefficient and temperature. In 
contrast a ferromagnetic coherence length which is also a 
depth of the condensate penetration into a ferromagnct 
in FS system is much shorter = ^ D /h provided the 
exchange energy h is rather large T which is usually 
fulfilled. 

Secondly, the penetration of Cooper pair wave function 
into the ferromagnetic region (F) is characterized by the 
damped oscillatory behaviour of a correlation function 
/ ~ exp(— x/^i?) sin(a;/^i?) which is a result of exchange 
splitting between energy bands of conduction electrons 
with different spin projections. In fact the origin of os- 
cillations is the same as for the Fulde- Ferrcl- Larkin - 
Ovchinnikov stated. This results in many new effects, 
such us spatial oscillations of the density of states^, a 
nonmonotonic^ or reentranljSi behaviour of the critical 
temperature as a function of a ferromagnetic layer thick- 
ness in layered FS structures. Also it is responsible for 
the formation of Joscphson tt junctions- and spin valves^. 

Despite the short coherence length in the ferromag- 
netic region there is a possibility of a long range proxim- 
ity effect in FS structures with inhomogcneous magnetic 
structure. In experiments on FS systems with strong 
ferromagnets an anomalously large increase of the con- 
ductance below the superconducting critical temperature 
Tc was observe d^i^i^° . Also recently the Andrcev interfer- 
ometer geometry was used to measure the phase sensitive 
conductance modulation in the FS system with helical 
magnetic structure^i. 

The first theoretical analysis of a long range proxim- 



ity effect in FS structure with inhomogcneous magne- 
tization was done for a Bloch-typc domain wall at the 
FS interfaced. It was shown that a superconducting 
correlation function contains components which survive 
at the distances of order of the normal metal correla- 
tion length from the superconducting boundary. These 
long range superconducting components have non-trivial 
structure in spin space. Conversely to the ordinary 
Cooper pairs which have a singlet spin structure they 
have a triplet spin structure which corresponds to corre- 
lations between electrons with the same spin projections. 
Therefore the long range superconducting components in 
FS systems arc usually called the long-range triplet com- 
ponents (LRTC). The LRTC can be generated in sys- 
tems with Blocbi^ and Neeli^i^^ domain walls or heli- 
cal magnetization patterni^. The long range proixim- 
ity effect was shown to exist in multilayered FS struc- 
tures with noncoUinear magnetization in different ferro- 
magnetic layersiSiiiii^ii^. Large attention has been paid 
to the investigation of long range Joscphson effect due 
to LRTC in FS systems [see Rcfl20l for a review]. Re- 
cently in Reffiol a multilayered SFIFIFS structure has 
been shown to demonstrate a controllable crossover be- 
tween long range triplet and short range singlet Joseph- 
son effects with the rotation of the magnetic moment of 
any of the F layers. 

The present paper is devoted to another possibility 
of controllable switching between long and short range 
proximity effects by employing the peculiar properties of 
ferromagnetic nanoparticles. In some sense the magne- 
tization of a nanoparticle is more simple than the do- 
main structure of macroscopic ferromagnets, therefore, 
theoretical findings could be proved by experiments with 
nanoparticles. It is now well-understood that a mag- 
netization distribution in a single particle is determined 
by the competition between the magnetostatic and ex- 
change energies. If a particle is small, it is uniformly 
magnetized and if its size is large enough a non-uniform 
(vortex) magnetization is more e nergy p referable (see, for 
example, Refs. HHIllIllIlilliiilaTliy) . Besides the ge- 
ometrical form and size, the state of the particle depends 
on many other factors. For example, if the ferromagnetic 
particle is initially in the vortex state then by applying 
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a homogeneous in-plane magnetic field one can shift the 
center of a magnetic vortex towards the particle edges^S. 
If the magnetic field is large enough the magnetic vor- 
tex annihilates, i.e. the particle becomes homogeneously 
magnetized. Conversely, applying magnetic field to the 
homogeneously magnetized particle in the direction op- 
posite to its magnetic moment one can force a nuclcation 
of magnetic vortex. Experimentally the shifting of mag- 
netic vortex is observed as a linear growth of the particle 
magnetic moment which saturates at the field of vortex 
annihilation. A transition from homogeneous to vortex 
state leads to a large jump of the magnetic moment so 
the magnetization curve of a ferromagnetic nanoparticle 
is in general highly hysteretio^. 

In practice superconducting correlations in a ferromag- 
netic nanoparticle can be induced in planar geometry by 
lateral superconducting junctions connected to the par- 
ticle. Obviously if the particle is homogcnously magne- 
tized then no long range correlations are induced and the 
proximity effect is short range. If the particle is in vor- 
tex state the situation is not so obvious and the special 
investigation is needed. Throughout this paper we will 
consider only the vortex state of the ferromagnetic par- 
ticle. We will show that for a circular particle there is 
no long range superconducting correlations if the mag- 
netic vortex is situated at the center. However if the 
magnetic vortex is shifted from the center by an exter- 
nal field there appear long range correlations. Moreover 
an axial anisotropy of geometric form of the particle also 
leads to a long range proximity effect. 

The structure of this paper is following. In the next 
section we describe our model, present the basic equa- 
tions and give a qualitative explanation of the long range 
proximity effect in a ferromagnetic nanoparticle with vor- 
tex magnetization. In Section III we present our main re- 
sults which are discussed in Section IV. Finally the con- 
clusions are given in Section V. 



II. MODEL AND BASIC EQUATIONS 

We consider a system shown schematically in Fig.([T|). 
It consists of a ferromagnetic nanoparticle and a lateral 
superconducting lead. The particle magnetization is as- 
sumed to form a magnetic vortex state. The structure 
of magnetic vortex is shown in Fig. (Hb). It can be 
described by the rigid vortex model proposed by Usov 
and Peschany^i and by Guslienko^S,. Within this model 
magnetization has a z component only inside the core 
region which size is determined by a ferromagnetic ex- 
change length Ip^c ■ Outside this region magnetization lies 
within xy plane. Typically the exchange length is quite 
small I ex ^ lOnm compared to the sizes of ferromagnetic 
nanoparticles R ~ lOOnm therefore we will neglect the 
vortex core region throughout this paper. Thus if the 
center of magnetic vortex is situated at the point r = a 
the magnetization distribution can be written in the fol- 




FIG. 1: Sketch of the system considered: (a) Ferromagnetic 
nanoparticle with vortex magnetization and attached super- 
conducting electrode, (b) Magnetic vortex at the center of the 
circular particle and the polar coordinate system, (c) Shifted 
magnetic vortex and the polar coordinate system with the 
origin at the vortex center. 

lowing form: 

M = g,„7\/o[np,zo], (1) 

where Zq is a unit vector along z axis, Wp = (r — a)/|r — 
a| and a is a shifting vector of magnetic vortex center 
with respect to the origin. In polar coordinate system 
(p, a) with the origin at the center of magnetic vortex 
(see Figd]:) the magnetization distribution ([T]) takes a 
simple form: 

M = (7mAfo(sina, — cosa). (2) 

The equation ([2|) describes the magnetization vector curl- 
ing around the center r = a in a clockwise (counterclock- 
wise) direction for — -fl(— 1). Further we will as- 
sume a clockwise direction of magnetization rotation (see 
FigHb)- Note that in case of a circular ferromagnetic par- 
ticle the shift of a magnetic vortex from the center can 
be directly related to the external magnetic field Hext as 
follows: 



where Xp is the ferromagnetic nanoparticle linear mag- 
netic susceptibility^^. The corresponding distribution of 
the effective exchange field acting on free electrons can 
be taken as h = ho'M./Mo, where ho is determined by the 
value of the exchange integral (see e.g. Refill). 

Our goal is to find a condensate Green function in 
the ferromagnetic particle induced by an attached super- 
conducting lead due to a proximity effect (see FiglT^). 
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We consider the "dirty limit" assuming that the mean 
free path of electrons is much shorter than all coherence 
lengths: I S,s,^n,^f- The most restrictive condition is 
I ^ since the ferromagnetic coherence length is much 
shorter than coherence lengths in superconductor and 
normal metal ^at. It imposes certain limitation on the 
magnitude of exchange interaction which means that the 
ferromagnetic should not be very " strong" . 

To analyze a proximity effect in ferromagnetic parti- 
cle we will use Usadel equations for quasiclassical Green 
functions. Following the scheme presented in detail in 
review^S we introduce a matrix Green function^ 
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Here g is normal and / is anomalous Green functions 
which are matrices in spin space. A space where the 
matrix g is defined is a Gor'kov-Nambu space. We will 
denote Pauli matrices in Gor'kov-Nambu space as and 
in spin space as (T^ (i = 1, 2, 3). Unit matrices are fo and 
(To correspondingly. Following Ref[23 the spinor basis for 
Green functions is taken in the following form: 



9 = 



/ 



TT Ti 

it a 



n TT 
a TT 



where t and [ denote the spinors corresponding to spin 
projections Sz = ±1/2. 

It is convenient to use a transformation of Green func- 
tion g suggested by Ivanov and Fomino\^ g = VgnewV'^ , 
where 



V = exp {j~i^{T3 - ^0)0-3) . 



(4) 



After this transformation is done the Usadel equation for 
the matrix Green function g takes the following form: 



where second term is small \F\ <^ 1. Then we obtain a 
linearized equation for F: 

DV^F - 2\uj\F - isgn{Lo){fQ(h ■a),F} = 0, (8) 

where {..} is anticommutator. The linearized boundary 
condition for the function F at the S/F interface is: 



n VF^ Fsh, 



(9) 



where 7 = Rb<y, while Rh is the interface resistance per 
unit area and a is the conductivity of the ferromagnet, 
n is a unit vector normal to boundary. The anomalous 
function in bulk superconductor is: 

Fs = (n sin tp- T2 cos (p) a^F^cs ■ 



Here F^cs = Aq/ ^Aq -f o;-^, where ip and Aq is the phase 
and module of the superconducting order parameter. 

Note that in Eq.([8]) the components of F proportional 
to fi and T2 are not coupled to each other. Thus in fer- 
romagnetic region anomalous function has the following 
structure: F ~ (fi sin(/? — f2 cos(y9) /. where / is a ma- 
trix in spin space. For the function / (when matrices in 
Namby space omitted) we obtain the following equation 
in ferromagnetic region: 

DV'^f - 2\uj\f - isgn{uj)\Y ■ {a, /} = 0. (10) 

The solution of Eg. pTH) can be found as a superposi- 
tion: 



/ = ao^o + aio'i + 020-2 + 030-3. 



(11) 



In this expansion the first term corresponds to the sin- 
glet component and the last three terms correspond to 
the triplet components of anomalous function with dif- 
ferent directions of Cooper pair spin. Note that after the 
transformation ^ the spin space basis for the anomalous 
function / can be symbolically written as follows: 



/ = 



TT 
TT 



-TT 
TT 



DV{gVg) - uo [f^, g] - i [f3(h • a),g] - [A, g] = 0, (5) Therefore it can be seen 



where [..] is a commutator, D is a diffusion coefficient, lo is 
Matsubara frequency and h ~ {hx, hy,hz) is an effective 
exchange field. The gap function is given by 

A = (fi/mA - f2i?eA)(7o, 

If there are no superconducting correlations in the nor- 
mal metal region then the Green function (in Matsubara 
representation) is given by 

g{uj) = sgn{uj)T3ao. (6) 

The Eq. ([5]) can be linearized assuming that 

g = sgn{uj)T3ao + F, (7) 



Si&i = 0, 

where Si is an operator of spin projection for a Cooper 
pair with respect to the i-th axis {i = x, y, z). If the vec- 
tor ftr — (01,02,03) is parallel to some real vector q in 
3Z? space then the Cooper pair spin projection on the vec- 
tor q is zero. This means that the Cooper pairs consist of 
electrons with the opposite spin projections, or in other 
worlds the spin lies in the plane perpendicular to vec- 
tor q. As we will see below the exchange field h collinear 
with the vector q effectively decouples the electrons lead- 
ing to the fast decay of Cooper pair wave function into 
the ferromagnetic region. Otherwise if the vector q (or 
more generally ftr) is not collinear to exchange field h 
the LRTC appear. 
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The equations for coefficients are: 

DV^ao - 2\uj\aa - isgn{uj)h ■ ftr = 0, (12) 

D\7^ai — 2|ti;|ai — isgn{uo)hxaQ = 0, (13) 

DV^a^ — 2\uj\a2 — isgn{uj)hyao = 0, (14) 



DV^oa — 2|w|a3 — isgn{uj)hzao ~ 0. 



(15) 



Now let us discuss tfie general structure of solutions 
of Eqs. m [HI [m [ISl)- If the magnetization and thus 
exchange field are homogeneous than it is easy to sec that 
there are two types of solutions of Eqs. (fT2l \T3\ [Ml fTS]) : 
(i) short range and (ii) long range modes. Indeed if the 
vector ftr = (ai, 02, 03) is parallel to the vector h then we 
obtain two equations for the functions oo and b = \ttr\- 

DV^aa-2\uj\ao-isgn{Lu)hb^O, (16) 



DV^b ~ 2\Lu\b - isgn{uj)hao = (17) 

which have solutions in the form: {ao,b) ^ exp(An • r), 
where A = ±(1 ± i)kh/V2 and kh = I/^f = \/h/D and 
n is a unit vector with arbitrary direction. These modes 
are short range ones since ferromagnetic exchange length 
is typically very short. One can see that short range 
modes consist of the singlet part of the anomalous func- 
tion with the amplitude given by coefficient gq. Also the 
is a nonzero contribution from triplet parts. The Cooper 
pair spin is directed perpendicular to the exchange field 
h. Therefore, such triplet superconducting correlations 
are suppressed by the exchange field on the same length 
scale as the singlet ones. 

On the other hand if the vector ftr is perpendicular to 
h then the Cooper pair spin can be oriented along h. In 
this case the destructive influence of exchange field on 
Cooper pairs is reduced. Indeed, from Eqs.p2l [T3l [T4l 
[T5|) we obtain that = and b satisfies the following 
equation 



D\/'^b-2\u;\b = 0, 



(18) 



which have a solution b ~ exp(An • r), where A = zbl/^jv 
and £,N = \/D/\uj\. These modes are long range ones 
because the coherence length in normal metal ^at can be 
rather large. Note that since Oq = these modes contain 
no singlet component, i.e. they contain only LRTC. 

In case of homogeneous magnetization long range 
modes can not be excited because of the zero boundary 
conditions for the triplet components: 



n • Vai = 



(19) 



for i = (1,2,3). The sources at the FS boundary exist 
only for a singlet component: 



where n is a unit vector normal to the boundary. How- 
ever it is not so for the inhomogeneous magnetization 
distribution. The well-known examples when LRTC can 
be excited are Bloch domain wall in a thin ferromagnetic 
wir&i^ or spiral magnetic structure which can be realized 
in some rare-earth metalsi^. Also recently the case of 
Neel domain walls in planar proximity FS structure was 
investigatedi^ii^. 

Now let us consider magnetic structure with large 
scale inhomogeneities. In zero order approximation for 
short range modes we obtain the Eas. (|16ll7p for uq and 
6 = |ftr| again, although the direction of vector a adaibat- 
ically depends on the coordinate: ftr = bh/h. The so- 
lution can be written in the following form: {ao,b) = 
{A, B)F(r), where A and B = A{kh/X)^ are constant 
and F{r) ~ exp(A • r). The boundary conditions p9l20p 
can be written as follows: 



n • Vao = Fbcs/l 
h(n- V6) = -b (n- V)h, 



(21) 
(22) 



where n is a unit vector normal to the boundary. There 
are two short range modes which decay far from FS 
boundary in the ferromagnetic region, say with Ai = 
kh{l + i)/V2 and A2 = kh{l - i)/V2. Taking the su- 
perposition of these modes with arbitrary coefficients Ai 
and A2 we obtain from Eq. ((22|) : 



Ai{Xi - Si) + A2iX2 - S2) = 

Ai{Xi - S2) + A2{X2 - Si) = 0, 

where 5*1 = {n-V)hx/h and S2 = {n-\/)hy/h. This linear 
system of the homogeneous equations has a solution if 
and only if Si = S2, i.e. 



(n • W)hx = (n • W)hy 



(23) 



n • Vao = Fbcsh, 



(20) 



This condition is fulfilled only in some special cases. The 
most trivial of them is a homogeneous magnetization dis- 
tribution. Another particular case when condition (j23p is 
fulfilled is that of a circular ferromagnetic particle if the 
magnetic vortex is situated at the center of the particle. 
Indeed in this case h^.y depend only on 9 and therefore 
(n • V)hx,y ~ {d/dr)hx.y = 0. Otherwise if the magnetic 
vortex is shifted from the center or the particle shape is 
axially symmetric the condition is not fulfilled. It 
means that taking into account the short range modes 
only one can not satisfy the boundary conditions and 
with necessity the long range modes are excited. 

The above qualitative description of the eigen mode 
structure is based on the assumption of adiabatically slow 
variation of magnetization and exchange field h. On the 
other hand in boundary condition (|22p appears a deriva- 
tive of h which in fact is a source for long range modes. 
Below we will find the corrections to the above adiabatic 
structure of short range modes. We will show that even 
if these corrections are taken into account it is still not 
possible to satisfy boundary conditions considering 
only the short range modes. 
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III. STRUCTURE OF SHORT- AND LONG- 
RANGE MODES IN MAGNETIC VORTEX. 

For further considerations it is convenient to introduce 
new functions b± = ai ± ia2 . Taking the magnetization 
distribution in the form ([2]) we obtain: 

(V^ - kl) ao-i'$sgniu;) [S*{r)b+ + S{r)b^] = 0, (24) 



(V^ - kl) b+ - isgn[u)klS[r)ao = 0, (25) 



isgn{uj)klS*{r)ao^Q, (26) 



where fc^ = 2\u)\/ D and k\ = ho/D. We have introduced 
the foUowing function: iS'(r) = [hx + ihy)/ho, where ho = 



hi + hi. 



A. Short range modes. 

Usually the ferromagnetic coherence length = 1/fc/i 
is very short. Most importantly it is much smaller than 
the size of a particle R and the characteristic scale of the 
magnetization distribution given by the function S(r). 
Therefore solutions of Eas. (|24l25l26p with effective wave- 
length can be described within quasiclassical approxi- 
mation. Also we neglect terms proportional to k^. Phys- 
ically it is justified since usually the normal metal coher- 
ence length ~ l/kui is much larger than the ferromag- 
netic coherence length £^p. 

Then we obtain the solution of Eqs. (|24l25l26p in the 
following form (see Appendix \X\ for details): 



a,o(p)e'^«, b+ -- 
Then we obtain: 



"1 d 




A'P 1 

_ h2 






P 



apa~sgn{uj)Y {bp- 



%+) = 
(30) 



0, 



dp 



ld_ 

pdp 

ld_ _ 

pdp \dp 



{M + If 



kl 



{M-lf 



bp+~sgn{uj)klapo = 0, 
(31) 

bp-+sgn{uj)klapo = 0. 

(32) 

The behaviour of solutions of Eas. (|30l31l32p depends 
on the ratio of the ferromagnetic particle size and nor- 
mal metal coherence length R/S,n- Indeed if i? ^ ^at 
these modes decay at the length ^jv- This is not very 
interesting case both for the experiment and for the the- 
oretical study. Another limit which can be investigated 
analytically is realized when ^ R. It means that 
the decay of the long range modes on the size of a fer- 
romagnetic particle is weak. This condition is the most 
favorable for investigation of the long range proximity ef- 
fect. Therefore we neglect terms proportional to fcj from 
Eas. (|30l3ll32l l. 

It is possible to find the long range modes as expansion 
by the orders of small parameter {Rkh)~^. The details 
of the calculations are shown in Appendix [Bl We obtain 
the following solution: cipo = 0, 



bp+=bp^=Bp^^\ 
where B is an arbitrary coefficient. 



(33) 



IV. RESULTS. 



ao = F{9) exp(A n • r) 



(27) 



1,2 / 2 

6+ = iF{0)sgn{uj) exp(A n • r) ^ fl - - (n • V) ) S*, 



(28) 



b- = iF{e)sgn{uj) exp(A n • r)-| ( 1 - -(n • V) ) 5 



(29) 

where F{d) is an arbitrary 27r periodic function and A = 
Ai,2 -fc,,(l±i)/y2. 



B. Long range modes. 

Now we are going to consider slow modes of 
Eas. (|24l25l26[) . For this purpose we choose the co- 
ordinate origin at the magnetic vortex center (p, a) 
(see FigHJ;). Then we have S'(r) = — ie'" and there- 
fore Eqs. (|24l25l26p allow separation of variables: = 



We will find the distribution of anomalous Gor'kov 
function in a ferromagnetic nanoparticle induced by a su- 
perconducting electrode which is attached to the particle 
as it is shown in Fig.([T]). The superconducting electrode 
attached at some point to the ferromagnetic sample can 
be modeled by the angle-dependent transparency of the 
FS interface 7 = ^^{6) in the boundary condition (|20p . 
For simplicity we can consider a Gauss form of trans- 
parency: 



7 = 7o exp 



{{e-Oo) mod 2TrY 



(34) 



where SO is determined by a junction width d = 
R{5e /2tt). 

Let us start with a general consideration. The bound- 
ary conditions for the coefficients Oq, &+, 6_ at the bound- 
ary of a ferromagnetic particle read: 



n • Voo 



Fu, 



l{0) 



(35) 
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n • V6± = 0. 



(36) 



To satisfy the boundary condition for ao we take the 
superposition of solutions (|27l28l29p corresponding to 
Ai = kh{l + i)/V^ and Ai — fc/i(l — with arbi- 

trary functions Fi,2(^)- We will take into account only 
those solutions which decay far from the FS boundary. 
Using the expression ^ for the vortex magnetization and 
taking into account that S{r) = — ie*" from Eqs. p5l36p 
we obtain: 



Fbc 



(37) 

hn- V&,+ = 
(38) 

-n- Vfe(_ =0, 
(39) 

where &z±(r) are the contributions of the long range 
modes. The structure of the long range modes yields the 
following relation for the coefficients e~'"6/_)_ ~ e^"bi_. 
Let us denote sgn{uj)hn^ — e^'"5/_|- ~ e*"6;_. Then from 
Eas. (|38l39p we obtain: 



1- — (n- Va) 
A," 



1 + — (n- Va) 

A, 



Fi ^ F2 
Ai A2 



n • Va 



70- 



Fi F2 



(n-Va)(^ + ^)=z 
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(40) 



(41) 



One has Ai,2 ~ fc/i, therefore the r.h.s. of Eq. ((40|) is 
small and with good accuracy Fi/Xi + F2/X2 = 0. The 
Eas. (|37l4ip then yield 



1,2 



Fte,/(27Ai,2) 



and 



n- V6, 



i F 



bcs 
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V2kh 7 



(n • Va) 



(42) 



(43) 



We search the contribution from the long range modes 
as a superposition: 



bio = E 



(44) 



where p, a are the polar coordinates relative to the cen- 
ter of a magnetic vortex. We use numerical methods to 
calculate the coefficients in the sum (|44p . We assume the 
angle dependent transparency in Eq. ([55)1 in the form ((M]) 
with S9 = 0.02 and the value of the ferromagnetic coher- 
ence length — Q.02R. Further we will consider two 
typical cases: (i) magnetic vortex in a circular particle 
shifted from the center of it and (ii) magnetic vortex at 
the center of a particle having the elliptical shape. 




FIG. 2: Amplitude of the triplet anomalous function in- 
duced due to the proximity effect. The position of mag- 
netic vortex center is marked by the white circle, (a) Shifted 
magnetic vortex in circular particle, vortex shifting vector is 
a = (0, —R/2); (b) magnetic vortex at the center of elliptical 
particle with axes ratio Rx/Ry = 1.5. 



A. Shifted magnetic vortex. 

Let us assume for simplicity that the shifting vector is 
directed along x axis: a = [a^, 0). The vector normal to 
the boundary is directed along the disk radius: n = r/r. 
Then the short range modes are given by Eas. (|27l28l29p 
with 



S{v) 



i{a — re^^) 



(r^ — 2arcos0 + a'^) 



2^1/2 



and 



(n- V)S' 



dS 
Ik 



(r ^ — 2ar cos 9 + a^) 



2^3/2- 
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The boundary condition for the long range modes ([43 
takes the foUowing form: 

dbia. _ i Fbc 
dr 1^=^ - V2fc. 7 

where 



-Q{0), 



(45) 



The boundary condition for the fong range modes (j43p 
then takes the form ((45|) with 



Rl-Rl i?2cos2 
Q{e) = sin(2^) " ^ ' ^ 



R-^y C0S2 ( 



i?,? sin^ ( 



Ri sin^ 



a sm ( 



(i?2 - 2ai?cos6i + a2)' 

In general, the amplitudes of the short range modes 
given by Eas. (|27l42p are determined by the dimension- 
less factor ^f/io- From the Eg. ([55]) it is easy to see 
that dbio/dr^r = i?) ~ ho/R. Thus, when the vortex 
shifting is small (a ^ R) the amplitude of LRTC is de- 
termined by the dimensionless factor (^p /jo){a/ R), i.e. 
it is [R/a) times smaller than the amplitudes of the short 
range triplet components. 

In case when a junction with a superconducting lead 
is narrow 56 ^ 2tt, the amplitude of the LRTC is deter- 
mined by the function |Q(0o)|- One can see that the max- 
imum amplitude is obtained when cos^o = 2aR/{R'^ + 
a^). On the other hand the long range proximity effect is 
absent if = or tt. This is caused by the symmetry of 
the magnetization distribution. In such case the magne- 
tization is constant along the direction of surface normal 
vector at the point where the superconducting lead is at- 
tached. Therefore there appear no source for LRTC at 
the FS boundary. 

To demonstrate the enhancement of the LRTC in the 
ferromagnetic particle with the shifted magnetic vor- 
tex we plot in Fig|2^ the distribution of the amplitude 
of the triplet part of the anomalous function \ftr\ = 
A/|ai|2 -I- |a2|2 [see expansion ([TT]) ]. We choose the posi- 
tion of a superconducting contact = and the mag- 
netic vortex shifting vector a = (0, ay). 



B. Magnetic vortex in elliptical particle 

Now let us consider the situation when the mag- 
netic vortex is situated at the center of a particle but 
the particle itself has elliptical shape. The boundary 
of the elliptical particle is determined by the equation 
{x/ Rx)^ + {y / Ry)^ = I. It is convenient to write the 
vector normal to the boundary in the polar coordinate 
frame n = rirfo + ne9o where 

Rl cos2 e + Rl sin^ 9 



i?4 cos2 e 



Ri sin^ e 



no 



{Rl ~ i?2)sin(26') 
2,/i?^ cos2 e + i?4 sin2 61 



Then the short range modes are given by Eqs. (j27l28l29p 



with S'(r) 



-le 



and (n • V)5 = {n0/r){dS/d9), or 



(n.V)5 = e^Mn(20)^-^. 

ZUxKy 



IRl cos^ e + Rl Bin^e 



i?4 cos2 ( 



One can see that when the shape of the particle is 
nearly circular 6R ■ 



'|i?2 - i?2| <^ Rx.Ry the amph 
tude of LRTC is determined by the dimensionless fac- 
tor {£_f/7o){SR/Ro), where SR = 



i?,2| is a mea- 
sure of axial anisotropy of the elliptical ferromagnetic 
nanoparticle and Rq = ^ R^ + Ry- 

Since the center of the magnetic vortex is assumed to 
coincide with the particle center we search the long range 
modes in the form of expansion (|44p with p = r and 
a = 9. Then we obtain: 



n • V6;c 



m? + 1 + ingm 



Going along the same lines as in the previous section 
we find the coefficients Cm numerically and obtain the 
distribution of the amplitude of the triplet component of 
the anomalous function shown in Figl^b- We choose the 
position of the superconducting contact 9q ~ i: /A. 



DISCUSSION 



,x 10 



o 




FIG. 3: A superconducting phase sensitive correction to the 
local conductance as a function of the magnetic vortex dis- 
placement with respect to the center of a circular ferromag- 
netic particle (see the insert). Different curves correspond to 
the angle x values (from bottom to top): X = 0' 



2n Btt 47r 
^' 10 ' 10 ' 10 ' 2 • 



Ri sin^ 9 



Experimental observation of the proximity effect in 
FS structures can be done for example using Andreev 
interferometer geometry to measure the modulation of 
the conductivity of a ferromagnetic sample as a func- 
tion of the phase difference between the superconducting 
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leadsii. Therefore it is interesting to investigate the in- 
fluence of the long range proximity effect on the trans- 
port properties of ferromagnetic nanoparticles. Let us 
consider a system shown in Figl3] (see the inset). We 
assume that there are two superconducting leads with 
different phases of the superconducting order parameter 
Lpi,2 attached at the different points to the circular fer- 
romagnetic nanoparticle. The normal lead measures the 
conductance of the system. In case of a point junction 
with normal lead one can use a general relation between 
a zero-bias tunneling conductance G and a local density 
of states (LDOS) v in the ferromagnetic particle at the 
junction point: 

G = G„(z/K), (46) 

where G„ and i/„ arc the point junction conductance and 
LDOS in the normal state of the ferromagnetic particle. 
The above expression for the local tunneling conductance 
is valid only if the voltage drops in the small vicinity of 
the junction point. This condition can be obtained as- 
suming, for example, that the potential surface barrier is 
so high that all voltage drops just at the interface between 
the normal lead and the ferromagnetic particle. But in 
case of a point junction Eq. ((l^ can be used even for an 
ideal interface because the voltage drops at the distance 
determined by the junction size. Note that it is not so 
if, for example, a conductance of one-dimensional wire is 
consideredi^i^i. We will assume that the junction size is 
much smaller than other characteristic lengths and cm- 
ploy the expression (|46p for the tunneling conductance. 

Having found the condensate function /, wc can calcu- 
late the LDOS in the ferromagnetic region. The LDOS 
is given by the general formula^^: 

where to = —ie + and the trace is taken in both the 
Gor'kov-Namby and spin spaces. Using the normaliza- 
tion condition cf^ + ff~^ = 1 and the smallness of the 
condensate function, we obtain the correction to the con- 
ductance of the point junction: 

<5G/G„ = -ii?eTr(//+). 

The anomalous function F has the following structure in 
Gor'kov-Nambu space: 

F ^^fi (sin ipiTi - cos (piT2) , 

i 

where Ai 2 = Ag exp(i(/3i 2) are the gap functions in the 
superconducting leads. Therefore 

and 



Thus we obtain: 

SG/Gn = ^ReTr (/^ + /| + 2A/2 cos ^) , (47) 

where ip = ipi — ip2. Probably the most important for ex- 
periments is the conductance correction in Eq. (j47p which 
depends on the phase difference ip due to the interference 
between the anomalous functions induced by different 
superconducting leads 5G/G„ = cos (pReTr{f if 2)- In 
FiglHlwe show the dependence of the amplitude of con- 
ductance modulation SG = G„i?eTr(/i/2) on the dis- 
tance of the magnetic vortex center from the center of 
the ferromagnetic particle. Different curves in this plot 
correspond to the different directions of vortex shifting 
vector a (see the sketch of the system considered on the 
insert in Figl2]). We normalize the conductance to the 
following value Go = G„(7o/Cf) which is entirely deter- 
mined by the fixed parameters of the system. 

Analyzing Fig [3] one can see that the strongest effect 
is achieved by shifting the vortex symmetrically with re- 
spect to the superconducting contacts (x = 0). On the 
contrary, the effect of conductance modulation is very 
small in case of the vortex shifting along the line be- 
tween two superconducting leads (x = '"'/S, top curve). 
As we have discussed above, in this case the LRTC are 
weak due to the symmetry of the magnetization distri- 
bution. A non-zero value of the conductance modula- 
tion in this case is caused only by the finite width of 
the superconducting junctions used in calculations. Fur- 
thermore, in FigO all the curves, except for the top one 
which corresponds to x = 7r/2, demonstrate strong asym- 
metry with respect to the sign of the vortex displacement 
SG{a) ^ SG{—a). Such asymmetry is caused by the sys- 
tem geometry, since we consider a conductance of only 
one point junction. As one can see if the magnetic vortex 
shifts towards the normal contact (positive a in Fig[3]) the 
conductance modulation appears to be very small com- 
pared to the case when the magnetic vortex shifts in the 
opposite direction (negative a in Fig[3]). This effect can 
be understood if we recall that the long range modes are 
strongly suppressed near the vortex center [see Eq. ([55]) ]. 
Thus even if the overall amplitude of LRTC is increased 
with |a|, the local value of anomalous function at the 
junction point is decreased if the magnetic vortex center 
shifts towards the junction point. 

The shift of the magnetic vortex is unambiguously de- 
termined by the magnetic field [see Eq. ([3])], therefore 
the asymmetry 5G{a) 7^ SG{—a) will be revealed in the 
conductance dependence on the external magnetic field: 
SG{H) 7^ 6G{—H). But in reality one always has two 
contacts and the total conductance correction is a sum 
of the contributions from each contact. Thus the result- 
ing behaviour of the conductance should depend on the 
position of the points where superconducting and normal 
contacts are connected to the ferromagnetic particle. In 
particular, if the system geometry is symmetric with re- 
spect to the spatial inversion the conductance correction 
will not depend on the sign of vortex shifting a as well 
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as on the sign of the magnetic field 5G{H) = 6G{—H). 

In Fig[3] the modulation of conductance is shown not 
for the entire range of the magnetic vortex displace- 
ments from the particle center. The reason is a growing 
complexity of numerical calculations because when the 
magnetic vortex center approaches close to the particle 
boundary one has to take into account too many angu- 
lar harmonics in the expansion (|44p . We expect further 
monotonic growth of |(5G'(a)| until \a\ < R. If the vortex 
displacement distance becomes larger than the particle 
radius \a\ > R, the vortex actually leaves the particle. 
Such magnetization state often is referred as "buckle"—. 
Further increase of |a| describes in fact a continuous tran- 
sition to the homogeneously magnetized state. Therefore, 
the conductance correction should eventually vanish as 
|a|^0. 

The overall magnitude of the conductance modula- 
tions is determined by many factors. One of them is 
a vortex displacement, which can be regulated by the 
external magnetic field. Other factors are determined 
by the geometry of the system, e.g. width of super- 
conducting leads and the particle size R. Also there 
is a dimensionless factor ^f/7oi which depends on the 
material parameters: ferromagnetic coherence length S^p 
and 70 = Rpo'int, where Rp is the resistance per unit 
area of PS interface and aint is the conductivity of 
ferromagnetio^S. This factor determines the amplitude 
of the anomalous function within the ferromagnetic re- 
gion and should be small within our calculation scheme, 
because we consider the linearized Usadel equation. For 
a particular configuration shown on the inset in FigOwe 
obtain the maximal amplitude of conductance modula- 
tion SG ~ 10~^(^f/7o)G'„, where Gn is the unperturbed 
conductance in the normal state of the particle. Taking 
for example ^f/io ^ 10~^ we obtain that 6G ~ 10~'*G'„. 
To have a better effect in experiment one should try to 
increase the ratio S,f /lo- For example this can be ob- 
tained by using not very strong ferromagnetic material 
with relatively large e.g., Cu-Ni alloys^, characterized 
by rather large coherence lengths: S^f ^ lOnm. However 
the magnetic vortex has been observed in rather strong 
ferromagnets such as Co or Pe with ^f ~ Inm. On the 
other hand, one can try to improve the properties of the 
superconducting contacts, i.e. to use the contacts with 
low interface resistance Rf- 



VI. CONCLUSION 

To summarize we have investigated the proximity ef- 
fect in the ferromagnetic nanoparticle with nonhomoge- 
neous vortex magnetization distribution. We have de- 
rived a general solution both for the short range com- 
ponents and the long range triplet components of the 
anomalous function. Quite generally it is shown that the 
long range proximity effect can be realized if the axial 
symmetry of the magnetization distribution is broken ei- 
ther due to the shifting of magnetic vortex with respect 



to the particle center or due to the angular anisotropy of 
the particle shape, which can be, for example, ellyptical 
in real experiments. Also we have considered the su- 
perconducting phase-periodic oscillations of the particle 
conductance in Andrecv interferometer geometry, which 
has been used recently to study the proximity effect in 
a conical fcrromagnet^i. We have shown that the am- 
plitude of conductance oscillations strongly depends on 
the direction of external magnetic field which determines 
the shift of magnetic vortex with respect to the particle 
center. For a particular case of a circular ferromagnetic 
particle the conductance oscillations are the largest when 
the vortex shifting is symmetric with respect to the su- 
perconducting contacts position. However, we suppose 
that the optimal direction of vortex shifting for the obser- 
vation of the long range proximity effect should depend 
on the system geometry, such as particle shape and posi- 
tion of the points where the superconducting and normal 
contacts are connected to it. 
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APPENDIX A: DERIVATION OF THE SHORT 
RANGE MODES 

We search quasiclassical solutions of Eqs. (|24l25r26p in 
the following form: 

(ao, &+,&-) = exp(A n • r)(ao5, 6.,), (Al) 

where n is a unit vector, A is large and functions 
aoq,b^q,b-q are slow. Note that in principle the direc- 
tion of vector n is arbitrary and should be determined 
from the boundary conditions. But we assume from the 
beginning that the spatial scale of the anomalous func- 
tion variation along the boundary is much larger than 
1/|A|. Thus we can consider vector n as a normal to 
the boundary of a ferromagnetic. Then at first we need 
to find A. Substituting functions in the form (|Aip into 
Eas. (|24l25l26|l we obtain: A* = -k^, i.e. A = {-ly^^kh 
which corresponds to the short range modes and A = 
which we will discuss. For quasiclassical envelopes we 
obtain the following equations: 

2A(nV)ao,+A2aoq-z^sgn(w) [S*{r)b+g + 5(r)6_,] ^ 0, 

(A2) 

2X{nV)b+q + X^b+g - isgn{uj)klS{r)aoq = 0, (A3) 
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2A(nV)6_5 + \%-q - isgn{uj)klS*{v)aoq = 0, (A4) 

Since all A have large real parts all the solutions decay 
or grow very fast. We will take into account only those 
which decay far from the boundary of the ferromagnetic 
particle. Then we should leave Ai = kh{\ + i)/^/2 and 
A2 = fc?i(l — Let us now find the solutions of qua- 

siclassical Eas. (|A2IA3IA4|) . We wiU use a perturbation 
method. 

Let us at first assume that a^q — const. Then to the 
zero order: 



1 d 



p dp \ dp 



d 



(M-1) 



21 



bp- + sgn{uj)klapo = 0. 



(B3) 

It is convenient to rearrange these equations introduc- 
ing new functions bs = bp^ + bp_ and b^i = bp^ — bp^ : 



^ + -^-^ja.o + ..nHf6. = 0, (B4) 



fc2 

= iaQqSgn{uj)^S{v) 



fc2 

6-g = iaoqsgn{uj)^S* [y) 



(A5) 



(A6) 



Note that we also can assume aoq = G(r). where G{y) 
is arbitrary but rather slow function. In this case two 
other coefficients b+q and b-q are proportional to G(r). 
The condition (nV)G' ^ |A| guarantees that this will 
not change the structure of eigcn modes. Substituting 
expressions (jA5IA6|) into Eas. (|A3IA4[) we obtain the first 
order perturbations 

b+q = -2iaoqSgn{u))j^{nV)S{r) 



^2 

-2iaoqSgn{u;)^{nV)S*ir). 



APPENDIX B: DERIVATION OF THE LONG 
RANGE MODES 



/ 92 



ydp^'-pd-p-^y^-l^''^ = '^ 

dp"^ p dp p^ 



bd :^bs-2sgn{uj)k,^apo = 0. 

P 

(B6) 

We will find the solutions of these equations as expansion 
by the orders of the small parameter {pkh)~^ assuming 
that the distance from vortex center is much larger than 
the ferromagnetic coherence length p^ £,f- It is easy to 
see that if fc/i — > 00, we obtain that 6^ = and Upo ~ 
and 
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1 d AP + 1 



dp"^ p dp 



0. 



(B7) 



The solution of this equation is bs — BgoP^^^^ ■ Then 
from the Eq. (jB6P we get: 



-sgn{uj) 



2M 



The long range modes can be found solving Eqs. 
(|30|31|32p with neglected terms proportional to A:^: 



B 



sO- 



ld_ f d_\ 
pdp\dp) p2 



I d f d 



pdp \ dp 



ipo - sgn{uj)-^ {bp- - bp+) = 0, 

(Bl) 



(A/ + 1)2 



bp+ - sgn{LL>)klapo = 0, 



(B2) 



The function bd is obtained from Eq. (|B4p : 
4M(5-4^/lPTT) ^^ 



Bso- 



Substituting it to the Eas. (jB5[) we obtain the next correc- 
tion for bg of the order [khp]^'^ which can be neglected. 
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